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XY and Heisenberg spins, subjected to strong random fields acting at few points in space with 
concentration Cr ^ 1, are studied numerically on 3d lattices containing over four million sites. 
Glassy behavior with strong dependence on initial conditions is found. Beginning with a random 
initial orientation of spins, the system evolves into ferromagnetic domains inversely proportional to 
Cr in size. The area of the hysteresis loop, m{H), scales as c^. These findings are explained by 
mapping the effect of strong dilute random field onto the effect of weak continuous random field. 

Our theory applies directly to ferromagnets with magnetic impurities, and is conceptually relevant 
to strongly pinned vortex lattices in superconductors and pinned charge density waves. 
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The random-field problem has a long history. In the 
early 1970s, Larkin argued [T] that random pinning of 
flux lattices in superconductors leads to a finite trans¬ 
lational correlation length which determines the criti¬ 
cal current |2]. A more general statement was made by 
Imry and Ma [3], who argued that a static random field, 
regardless of strength, breaks the system into domains 
below d = 4 spatial dimensions. Besides superconduc¬ 
tors, the concept of Imry-Ma domains was applied to 
random magnets [MS], disordered antiferromagnets |7|, 
spin-glasses [5|, arrays of magnetic bubbles [Ij, charge- 
density waves (CDW) [TnHT3] . liquid crystals |T3], the 
superconductor-insulator transition m, and superfluid 
3He-A in aerogels na HZ!- In the early 1980s the va¬ 
lidity of the Larkin-Imry-Ma argument was questioned 
by the renormalization group treatment USKH]. Scaling 
and replica-symmetry breaking methods [10HS0|> well 
as the variational approach Eiiisa, suggested a power- 
law decay of correlations at large distances, indicating a 
type of topological order (so called “Bragg” or “elastic” 
glass) which is robust against weak collective pinning. 
It was argued [33] that the energy associated with vor¬ 
tex loops prevented the XY Sd random-field system from 
completely disordering. The interpretation of decoration 
[31 [55] and neutron diffraction m experiments on flux 
lattices has been hampered by the fact that for weak 
disorder the correlation length can be very large, mak¬ 
ing it difficult to distinguish between large defect-free do¬ 
mains and a Bragg glass. Numerical work on Id [37] [3H]), 
2d [33] [10]), and 3d [SHU systems with quenched ran¬ 
domness established strong non-equilibrium effects, such 
as hysteresis and dependence on initial conditions. Re¬ 
cently, it has been argued and confirmed numerically 
on lattices containing up to one billion sites [11 that 
long-range correlations in the presence of a weak random 
field are controlled by topology. For a fixed-length n- 
component order parameter in d dimensions the behavior 
of the system is fully reversible, characterized by the ex¬ 
ponential decay of correlations, at n > d-|-1. Meanwhile, 


at n < d -I-1 topological defects leads to metastabillity 
and glassy behavior that depends on the initial condition. 

In general, for spin systems it is hard to come up 
with a physical mechanism that would generate a ran¬ 
dom field (different from magnetic anisotropy) linearly 
interacting with a spin at every lattice site. An exception 
is anisotropic antiferromagnet with fluctuating exchange 
interaction in a uniform magnetic field [7]. However, if 
disorder due to, e.g., magnetic impurities is confined to 
few points in space, then the exchange interaction — J^-S-s 
of the lattice spin s with the nearby impurity spin S 
is equivalent to the interaction with the magnetic field 
h = JyS. The model of a ferromagnet with dilute impu¬ 
rity spins with frozen random orientations is, therefore, 
equivalent to the random-field model in which the field 
only acts on spins located at few points in space. The XY 
spin model of this kind was recently studied by Okamoto 
and Millis m in connection with the problem of strong 
pinning of the CDW by dilute pinning centers. The com¬ 
puted correlations extended far beyond the average dis¬ 
tance between strong pinning centers. This study was 
inspired by the experimental evidence of extended phase 
correlations in the CDW in NbSe 2 m that disagreed 
with earlier theoretical predictions [IDHIS]. 

Our study of XY and Heisenberg spins on 3d lattices 
containing over four million sites is focused on metasta¬ 
bility. The measure of metastability is the area of the hys¬ 
teresis loop. We find that it scales as c^, where c^. <C 1 is 
the concentration of spins subjected to strong dilute ran¬ 
dom fields. If one starts with random initial conditions 
(RIC) for the spins the system evolves towards short- 
range order with exponential decay of correlations and 
with the correlation length inversely proportional to c^. 
In accordance with the findings of Okamoto and Millis, 
as well as with experiments on NbSe 2 , the correlation 
length that we obtained is much greater than the aver- 
age distance, 1/cr' , between the impurities. We provide 
a simple physical explanation to these findings. 
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XV model, Cr = 0.09 



Figure 1. Color online: Snapshot of XY spins in one layer of 
a 3d lattice after relaxation from random initial orientation in 
the presence of dilute random field. Spins whose orientations 
are frozen by the random field are shown in red. 

We consider the spin Hamiltonian 

7^ — 2 ^ ^ Jij^i ’ Sj ^ ^ hj • H • ^ ^ Sj, (1) 

ij i,r i 

with the nearest-neighbor exchange interaction, strong 
dilute random field h acting at r-sites with concentra¬ 
tion Cr 1, and the external field H, on lattices of 
a few millions spins of length s. In numerical work 
we consider infinite h, which implies that some spins 
with concentration Cr are frozen in random directions. 
Our numerical method combines sequential rotations of 
spins towards the direction of the local effective field, 
+ with energy-conserving 

spin flips: 2(sj applied with 

probabilities a and 1 — a respectively; a playing the role 
of the relaxation constant. This method is physically 
equivalent to slow cooling. Its fast convergence to the 
final state that is no longer relaxing has been demon¬ 
strated in Refs. 

We begin with evolution of the system from random 
initial orientation of spins. The randomly chosen Cr frac¬ 
tion of spins remain frozen in random directions, while 
other spins are allowed to rotate and relax to some fi¬ 
nal state in which the total magnetization is no longer 
changing. A snapshot of such a state is shown in Fig. 
Presence of topological defects and short-range order is 
apparent. The corresponding spin-spin correlation func¬ 
tion and its fit by the exponential are shown in Fig. 
This fit allows one to extract the ferromagnetic corre¬ 
lation length Rf. Its dependence on the concentration 
of impurities, Cr, is shown in Fig. The I/cr depen- 


Correlation Functions for XY model RIC, L = 160 
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Figure 2. Color online: Spin-spin correlation function (dashed 
lines) after relaxation from random initial orientation of spins 
in the presence of dilute random field. Solid lines provide the 
corresponding fit by the exponential. Distances are given in 
lattice units, with L being the size of the L x L x L system, 
a) 3d XY spin model, b) 3d Heisenberg spin model. 

deuce of Rf provides a good fit to the numerical data for 
both the XY model with two spin components and the 
Heisenberg model with three spin components. 

The state obtained by the relaxation from collinear 
initial conditions (with all spins initially oriented in one 
direction) always has non-zero magnetization. This is an 
indication of metastability and glassy behavior caused by 
the dilute random field. There may be as little profit in 
looking for the “needle-in-the-haystack” ground state of 
such a system as in looking for the ground-state config¬ 
uration of domains of a conventional permanent magnet 
with a macroscopic number of defects that pin domain 
walls. We, therefore, turn to the magnetic hysteresis as 
a measure of metastability. We compute the magnetiza¬ 
tion per spin ruz induced by the field applied in the Z- 
direction. Hysteresis loops for XY and Heisenberg mod¬ 
els are shown in Figs. [^and|^ The scaling in terms of 
H/c^ is very good, implying that the coercive field and 
the area of the hysteresis loop roughly scale as c^. 

Our findings can be explained in a simple manner 
by employing the Larkin-Imry-Ma argument. It is con¬ 
venient to rescale the problem to subvolumes of size 
Th ~ which represents the average distance be- 
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Rf vs Cr for XY model RIC, L — 160 



Rf vs Cr for H. model, RIC, L = 160 



Cr 

b 

Figure 3. Dependence of the ferromagnetic correlation length 
on the concentration of random field sites for the state ob¬ 
tained by evolution from random initial orientation of spins, 
a) 3d XY spin model, b) 3d Heisenberg spin model. 

tween the random-field sites. In the lattice units, the 
effective random field per spin inside such a volume is of 
order hr^, ^ h/r\. If we consider a group of spins ferro- 
magnetically correlated within a volume of size Rf Vh, 
the exchange energy per spin inside this volume is of 
order J/R'j, while statistical fluctuation of the Zeeman 
energy per spin is of order —hrf^jiRf. Mini¬ 
mization of the total energy with respect to Rf gives 
Rf ~ rf^(J/h)^ oc 1/cr, in accordance with our numerical 
results. In a similar manner one can compute the scaling 
of the coercive field. It roughly equals the field needed to 
overcome the effective field hRj, ~ hr^ /{Rf/rhY^'^ which 
aligns the spins on the scale Rf. Substituting and 
Rf from above, one obtains oc 1/r® ^ c^. This scal¬ 
ing of the coercive field and, consequently, of the area of 
the hysteresis loop is in agreement with our numerical 
results. 

Our model directly applies to a ferromagnet with mag¬ 
netic impurities having strong single-ion anisotropy that 
freezes their spins in random directions. We assume that 
such impurities interact strongly via ferromagnetic or 
antiferromagnetic exchange with the surrounding spins. 
Note that every magnetic system also has relativistic in¬ 
teractions, such as dipole-dipole and magneto-crystalline 


Hysteresis, XY model L = 160 



Scaled Hysteresis, XY model, L = 160 



Figure 4. Color online: Hysteresis loops for the 3d XY model 
with dilute strong random field for different concentrations of 
the random-field sites, Cr- a) Unsealed per-spin magnetization 
ruz vs H. b) Scaled per-spin magnetization ruz vs H/c^. 


anisotropy terms in the Hamiltonian, that would break 
the system into domains even in the absence of the spin 
impurities generating strong random fields. However, for 
sufficiently large concentration of magnetic impurities the 
formation of the short-range order would be dominated 
by their strong exchange interaction with the surround¬ 
ing spins rather than by the weak relativistic interac¬ 
tions in the system. In general, magnetic impurities must 
be of practical significance if the short-range-order ferro¬ 
magnetic correlation length, i?/ oc 1/cr, that they are 
responsible for, is smaller than the would-be average do¬ 
main size without the magnetic impurities. 

While the spin problem does not map exactly onto the 
flux-lattice problem or the CDW problem (see discussion 
of these issues in, e.g.. Refs. im and IT5]l , all three prob¬ 
lems are conceptually similar. Consequently, our study 
supports statements of Refs. HU and US that the cor¬ 
relation length in the CDW problem with diluted strong 
pinning centers is large compared to the average distance 
between the centers. Our prediction for the correlation 
length is Rf ^ Ijcr ^ Vh ^ t/cr ■ It would be interest¬ 
ing to test this prediction in NbSe 2 . It can also be tested 
on flux lattices that are pinned by strong pinning centers 
with concentration <SC 1. Scaling of the critical current 
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Hysteresis, Heisenberg model, L = 160 
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Figure 5. Color online: Hysteresis loops for the 3d Heisenberg 
model with dilute strong random field for different concentra¬ 
tion of the random-field sites, Cr- a) Unsealed per-spin mag¬ 
netization rriz vs H. b) Scaled per-spin magnetization ruz vs 

Hid. 


as jc'^ ^/R^f [2] would then imply jc oc d- 
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